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qq ■ We investigate continuous variable (CV) quantum teleportation using relevant classes of non- 

' Gaussian states of the radiation field as entangled resources. First, we introduce the class two- mode 

, squeezed symmetric superposition of Fock states, including finite truncations of twin-beam Gaus- 

£SJ ■ sian states as special realizations. These states depend on a set of free independent parameters that 

can be adjusted for the optimization of teleportation protocols, with an enhancement of the success 
probability of teleportation both for coherent and Fock input states. We show that the optimization 
procedure reduces the entangled resources to truncated twin beam states, which thus represents an 
optimal class of non-Gaussian resources for quantum teleportation. We then introduce a further 
class of two-mode non-Gaussian entangled resources, in the form of squeezed cat-like states. We 
O^l . analyze the performance and the properties of such states when optimized for (CV) teleportation, 

and compare them to the optimized squeezed Bell-like states introduced in a previous work [l2| . We 
discuss how optimal resources for teleportation are characterized by a suitable balance of entangle- 
ment content and squeezed vacuum affinity. We finally investigate the effects of thermal noise on 
the efficiency of quantum teleportation. To this aim, a convenient framework is to describe noisy 
entangled resources as linear superpositions of non-Gaussian state and thermal states. Although the 
presence of the thermal component strongly reduces the teleportation fidelity, noisy non-Gaussian 
states remain preferred resources when compared to noisy twin-beam Gaussian states. 

PACS numbers: 



c3 



' In recent years, much effort has been dedicated to the theoretical investigation and experimental production of 
ON \ highly nonclassical, non-Gaussian states of the radiation field Such states are endowed with enhanced nonclassical 
■ properties, such as entanglement and negative quasi-probability phase-space distribution, that can be useful for 
improving the efficiency of quantum information and communication protocols [l|, [H, 0, U, [H, . Importantly, a general, 
rigorous result states that, at fixed covariance matrix, some of these desirable properties are minimized by Gaussian 
states, which can then be considered as extremal "benchmark" or lowest-threshold states for quantum information 
applications Q. Therefore, at least in principle, the exploitation of non-Gaussian states and non-Gaussianity allows to 
quantitatively improve the performance of several quantum information protocols. For instance, non-Gaussian doners 
of coherent states turn out to be optimal with respect to Gaussian ones, exhibiting an enhanced single-clone fidelity [f| . 
Concerning quantum teleportation with continuous variables (CV), the success probability of CV teleportation can 
be greatly increased by using entangled non-Gaussian resources, " degaussified" through conditional measurements, 
consisting of photon-subtraction 0, 0, [HI, H2]- In particular, as shown in Ref. [l2|, the implementation of 
simultaneous phase matched multiphoton processes and conditional measurements can be used to introduce a general 
class of two-mode non-Gaussian entangled states, in the form of squeezed Bell-like states depending on an extra 
free, experimentally adjustable parameter. The optimization on this additional free parameter allows a remarkable 
increase of the teleportation fidelity for several relevant input states, compared both to Gaussian and other degaussified 
resources [l2j . From the point of view of quantum state engineering, several schemes for the generation of non-Gaussian 
states have been proposed 0,0, [01 Hi Hi Eli an d some of them, based on photon addition or subtraction, have been 
experimentally implemented to obtain photon-added and photon-subtracted single-mode states from Gaussian inputs 
fl&[l^ .l2fjl|. More recently, very advanced experimental realizations have been reported [2l|, HH, [23| . First, enhancement 
of entanglement and production of states with negative two-mode Wigner function have been successfully obtained 
via photon subtraction of a single delocalized photon from an initial two- mode Gaussian state [2l| . Remarkably, the 
photon-addition/subtraction operations, performed on a thermal light field, have also led to the demonstration of the 
commutation relation rules, one of the constitutive relations of quantum mechanics [22| . Moreover, using homodyne 
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detection and Fock-state resources, the experimental generation of large squeezed Schrodinger cat states has been 
achieved [23|. These optical cat-like states are of particular importance because they are strongly resilient against 
decoherence [24]]. 

In the present work, we aim to extend the results obtained in Ref. [HJ by generalizing the class of squeezed Bell-like 
states and by proposing an alternative class of non-Gaussian superposition states, and by showing that these two 
classes of entangled non-Gaussian resources allow for the optimization of CV quantum teleportation. Specifically, 
we investigate CV teleportation with two different classes of two-mode entangled resources: Squeezed symmetric 
superpositions of Fock states, and squeezed cat-like states. These two classes share two common features: Both 
arc associated to the squeezing operator as CV generator of quantum correlations, and for a trivial choice of the 
parameters, both reduce to the Gaussian twin beam (TwB). The last part of the work will be dedicated to the 
investigation of the effects of the presence of thermal noise on the performance of two-mode non-Gaussian states 
used as resources for quantum teleportation. In particular, we will consider the non-Gaussian resources obtained by 
superimposing squeezed Bell-like states or squeezed cat-like states and two-mode thermal states. Non-Gaussian states 
and thermal states are assumed to be independently generated, and overlap in a common spatial volume [25| . Due to 
the thermal contribution, the ensuing mixed states and their correlation properties are modified by the presence of 
thermal photons. We limit the discussion to the situation of ideal teleportation protocol, i.e. ideal Bell measurements 
and decoherence-free propagation in communication channels. Detailed analysis in the instance of the most general 
realistic situation, including various sources of noise, will be discussed elsewhere. 

The paper is organized as follows. In Scction[H we briefly recall the mathematical machinery needed for the description 
of the ideal Braunstcin-Kimble CV teleportation protocol in the characteristic function formalism. Next, we introduce 
the class of squeezed symmetric superposition of Fock states, and apply to them the optimization procedure defined in 
Ref. [l2l ]. In Section |TT1 we introduce the class od squeezed cat-like states, we apply the optimization procedure, and 
we compare results with those obtained in the previous Section. In Section [IIII wc will analyze the effect of noise on 
the performance of non-Gaussian resources for CV teleportation. Finally, in Section [IVI wc will draw our conclusions 
and discuss some perspectives for future research. 



I. CV TELEPORTATION WITH SQUEEZED SYMMETRIC SUPERPOSITION OF FOCK STATES 

The standard Braunstein-Kimble CV teleportation protocol, usually described in the domain of quantu m op tics by 
the Wigner function formalism [2(|, can be expressed equivalently in different mathematical approaches p7l.l28l.[29l.[30j 
( For a comprehensive review on CV quantum teleportation and information processing, see Ref. [3l|). Following a 
recent proposal, a very convenient description of the protocol can be given in terms of the characteristic functions 
of the various quantum states involved, including input, resource, and teleported states (30|. Denoting by pi n and 
Xini&in), respectively, the single-mode input state to be teleported and the associated characteristic function, and by 
P12 and Xi2(o!i , 0:2)1 respectively, the entangled two- mode resource and its characteristic function, it can be shown 
that the characteristic function Xout(ct2) of the teleported state takes the factorized form [30]: 

Xout{a-%) = Xin{oL 2 ) ^12(0:2 ,02) • (1) 

A commonly used measure of the success probability of a teleportation protocol is the fidelity of teleportation, i.e. 
T = Tr[pi n p ou t\. In the characteristic-function formalism, the fidelity reads 

T = i / d 2 \ X inWXout(-\)- (2) 

In the following we will adopt Eq. to analyze the efficiency of the CV teleportation protocol. 

In Ref. [IH, the squeezed Bell- like states have been exploited as non-Gaussian entangled resources optimizing the 

teleportation fidelity for various sets of input states. The general expression of squeezed Bell-like states reads: 

\4>)sb = [c? + C 2]- 1 / 2 5 12 (C){c 1 |0,0) 12 + e t %|l,l) 12 }, (3) 

where S 12(C) — e -c ' a i 2 + ^ aia2 is the two-mode squeezing operator, C = re 1 ^, \m\ ,m 2 )i2 = \ m i)i <8 |tt7. 2 ) 2 is a two- 
mode Fock state, and the c^s are real constants. A convenient paramctrization is c\ = cos 5 and c 2 = sin<5. At fixed 
squeezing parameter r, and at fixed phases <f> = 7T and 9 = 0, the entangled resource Eq. ([3]) is endowed with an 
additional phase S, that can be freely adjusted and thus used to maximize the fidelity Tsb'- 



F op t(r) = max Tsb(t,8) 





(4) 
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Such an optimization yields a remarkable enhancement in the success probability of teleportation for various input 
states [l2|. The exact analytical expression TsB{ r , 5) for the fidelity of teleportation of input coherent states using 
optimized squeezed Bell-like resources, i.e. Eq. (|Alj) . is reported in Appendix lAl The optimal phase S = 5 max , 
maximizing the fidelity Tsb{i", 5) is given by: 

8 ma x = \ arctan(l + er 2r ) . (5) 



It is now important to observe that the Bell-like contribution in Eq. ([3]) can be obtained from the truncation at first 
order (n = 1) of twin-beam states, i.e. <Si2(£)|0, 0)12 oc X)^Lo( — e%<P tanhs)"|n, n)i2, with £ = e lip s. This fact can be 
more clearly understood if one chooses, in Eq. ([3]), the paramctrization c± = 1 and C2 = tanhs (with ip = 9 + 71"). 
Thus, although being a non-Gaussian state, the squeezed Bell-like state can be extracted by a Gaussian squeezed 
state by means of a truncation procedure; in other words, it is a squeezed truncated-TwB. However, the connection of 
entangled non-Gaussian resources with the Gaussian twin beams may break down if we further generalize the squeezed 
Bell-like state © to the following class of non-Gaussian states, the two-mode squeezed symmetric superposition of 
Fock states: 

mssF = [c\ +4 + 4}~^ 2 S 12 (0{ci\Q, 0)12 + e^c 2 \l, 1) 12 + e ie *c 3 \2, 2) 12 }, (6) 

where Cj and 9i are real constants and phases, respectively. With respect to state (O, state ([5]) is not constrained to 
belong to the class of squeezed truncated- TwB. Namely, thanks to the additional four-photon term in the superposi- 
tion, the class of squeezed symmetric superposition of Fock states includes the squeezed truncated-TwB as a special 
realization. Now, it is very interesting to determine the subset of states \ip) ssf corresponding to an optimal fidelity of 
teleportation. Obviously, it is expected that the exploitation of the resource ^ will also allow a further optimization 
of the teleportation protocol. We compute the teleportation fidelity for an input coherent state \(3) and an input 
single-photon Fock state |1). A convenient parametrization of the CjS in Eq. ^ is provided by the hyperspherical 
coordinates in three dimensions: c\ = cos<5i, c 2 = sin<5i cosS 2 , C3 = sin^i sin^- The two-mode characteristic function 
XSSF of the entangled resource © can be calculated as 

XSSF(ai, a 2 ) = Tr[pDi(ai)D 2 (a 2 )] = sSF(i>\Di(ai)D 2 (a 2 )\ip)sSF , (7) 

where Di(cti) is the Glauber operator for the mode i. The explicit expression for Eq. ([7]) can be easily written by 
making use of the two-mode Bogoliubov transformations 

S\ 2 (C) a, S12 (C) = cosh r a, - e* sinh r a], (i # j = 1, 2) , (8) 

and of the relation 

(m\D(a)\n) = a m ~ n e~^ 2 L^(\a\ 2 ) , (9) 

where (•) denotes the associated Laguerre polynomial. The characteristic functions xc and Xf corresponding 

to an input coherent state \(3) and to an input single-photon Fock state |1), respectively, are given by: 

Xc(a in ) = e -h\-in\ 2 +2ii m[ain n } (10) 

XF(a m ) = e'^ a ^ 2 (1- |a m | 2 ) . (11) 

Given Eqs. (O, (TIT))) , and (fTl) . and using Eq. ([T]), the fidelities, defined by Eq. for the two different input states 
can be analytically computed. The resulting expressions are functions of the parameters r, <f>, Si, 61, 8 2l 6 2l but we 
do not report them, as they are too long and cumbersome. At fixed squeezing r, with maximally fixed phase <f> = n, 
the optimized fidelity is defined as 

Topt = max T{t,-k, Si, Qi,S 2 ,Q%) (12) 

where V denotes the set of independent free parameters V = {Si,8i,S 2 ,8 2 }. Thus, once the class of single-mode input 
states is specified, the optimization procedure of the available free parameters, in order to maximize the fidelity, leads 
to a sculpturing of the teleported output state. We compute the optimal fidelities (|12[) . for the above two classes of 
input states, by numerical maximization of the analytical expressions of T . First, we notice that the phases 6q and 9 2 
can be set to zero as they do not contribute to an improvement of T . The procedure leads to the following remarkable 
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result: Optimization over the set V yields maximal parameters compatible with a squeezed truncated- TwB as optimal 
entangled non- Gaussian resource. Thus, the optimal non-Gaussian resources for CV teleportation form a subset of 
the class dU) of the form: 



SSF 



[1 + tanh 2 s + tanh 4 s}~ 1/2 S 12 {-r){\0, 0) 12 + tanhs|l, l) 12 + tanh 2 s|2, 2) 12 }, 



(13) 



where, correspondingly to Eq. (|12[) . the real parameter s is fixed to the maximal value s = s, through the numerical 
maximization of the fidelity, i.e. T opt = max s T{r 1 s). In Fig. [1] we show the behavior of T ovt as a function of r for 
the input states \(3) (Panel I) and |1) (Panel II), with a squeezed symmetric superposition of Fock states as entangled 
resource. For comparison, also the curves T pti corresponding to the teleportation with a squeezed state and with a 
squeezed Bell-like state, are plotted. We observe that, as expected, a further enhancement of the fidelity is obtained, 
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FIG. 1: Optimal fidelity of teleportation J- op t, as a function of the squeezing parameter r, with cj> — n. In Panel I, we report the 
fidelity for the teleportation of input coherent states \j3) using different two-mode entangled resources: (a) Gaussian squeezed 
state resource (full line); (6) squeezed Bell- like state (dashed line); (c) squeezed symmetric superposition of Fock states (long- 
dashed line). In Panel II, we report the fidelity for the teleportation of input single-photon Fock states |1) using different 
two-mode entangled resources: (a) Gaussian squeezed state resource (full line); (6) squeezed Bell-like state (dot-dashed line); 
(c) squeezed symmetric superposition of Fock states (double-dot-dashed line). In plot I the value of (3 is arbitrary. The insets 
in both Panels give the maximal values of the parameter s = s as a function of r for given entangled resource and fixed input 
state. The plotstyles are chosen as specified above. 

compared to the case discussed in Ref. [l2j . 

In Ref. in order to understand the properties of the optimized resources, three quantities have been investigated: 
the von Neumann entropy E v n, in order to quantify the amount of entanglement present in the resource; the non- 
Gaussianity d n Q, in order to provide a measure of the non-Gaussian character of the resource; and the squeezed- vacuum 
affinity Q, in order to determine the degree of resemblance of the resource to a two- mode squeezed vacuum. In Fig. [21 
we plot E v n for the optimized non-Gaussian entangled resources Eqs. ^ and f) 13ft . and for the twin beam. It can 
be observed that all the curves exhibit very similar behaviors. At fixed r, and for a given input state, the optimized 
squeezed symmetric superposition of Fock states is the most entangled state, and the teleportation of a single-photon 
requires more entanglement than the teleportation of a coherent state. However, we notice that the behavior of the 
von Neumann entropies is in agreement with the behavior of the optimal fidelities; in fact, at fixed input state, a 
non-Gaussian resource with a higher teleportation fidelity is associated with a higher amount of entanglement content 
for any r.. Next, we examine the behavior of the measure of non-Gaussianity g?„g, which is defined as the distance 
between the considered state p and a reference Gaussian state pa, according to the following relation: 

Tr[p 2 ] + Tr[pj\ - 2Tr[pp G ] 
dnG = 2Trp] ' (M) 

where the Gaussian state pa is completely determined by fixing for it the same covariance matrix and the same first 
order statistical moments of the quadrature operators associated to state p [32| . In Fig. [3] panel I, we plot d n c for the 
entangled resources Eqs. © and fT3")) , optimized for teleportation. Looking at the behaviors of d n c corresponding to 
teleportation of the same input states, we see that, both in the instance of input coherent states and in the instance 
of input Fock states, the curves cross at a certain value of r. Therefore, according to the particular adopted measure 
of non-Gaussianity (|14p . a non-Gaussian resource with a higher teleportation fidelity does not exhibit, in general, 
a higher non-Gaussian character for any r. In order to quantify how a given state resembles a squeezed vacuum, a 



FIG. 2: Plot of the von Neumann entropy E v n as a function of r, with s = s (see the insets in Fig. [TJ, corresponding to 
the following entangled resources: squeezed symmetric superposition of Fock states, squeezed Bell-like states, and two-mode 
squeezed states, optimized for the teleportation of input coherent states \/3), and input Fock states The plot styles are 
chosen as specified in Fig. [I] 
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FIG. 3: Plot of the non-Gaussianity d n c (panel I) and of the squeezed-vacuum affinity Q (panel II) as a function of r, with 
s = s (see the insets in Fig. [1} , corresponding to the following entangled resources: squeezed symmetric superposition of Fock 
states, and squeezed Bell- like states, optimized for the teleportation of input coherent states |/3), and input Fock states |1). 
The plotstyles are chosen as specified in Fig. [1] 



measure of the squeezed-vacuum affinity can be introduced [12| . Such a quantity is defined as the following maximized 
overlap: 



g = max | 12 (-Clares (^)) 12 | 2 , 



(15) 



where | — £)i2 is a two-mode squeezed vacuum with real squeezing parameter — £, and |7/ves(f))i2 is any entangled 
two-mode resource that depends uniquely on the squeezing r as the only free parameter. In Fig. [3] panel II, we plot 
Q for the entangled resources Eqs. © and (fT5|) , optimized for teleportation. At fixed input states, the optimized 
squeezed symmetric superposition of Fock states exhibits a greater squeezed-vacuum affinity for any r with respect 
to the squeezed Bell-like state. This result is in complete agreement with the behaviors of the corresponding optimal 
teleportation fidelities and von Neumann entropies. Therefore, comparing the behaviors associated with resources 
optimized for the efficient teleportation of the same input states, we observe that the following hierarchy can be estab- 
lished: the enhancement of the teleportation fidelity corresponds to an enhancement of the entanglement content, and 
to an enhancement of squeezed-vacuum affinity. This hierarchy is not respected by the measure of non-Gaussianity 
d n G- This result shows that such a measure must be used and interpreted with care, case by case. As expected, in 
the plots of Fig. [3] panel I and panel II, for sufficiently large r, the curves, associated with the same resource, go to 
the same asymptotic saturation value. 

The following conclusions can be drawn from the above discussion: The squeezed truncated twin beam, with trun- 
cation at 77, = 2 (four-photon term), is an optimal entangled non-Gaussian resource suitable for sculptured quantum 
teleportation. We argue that this result can be generalized to any order n > 2 in expansions of the squeezing operator 
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at finite order. We thus conjecture that, in general, efficient entangled non-Gaussian resources for CV quantum tele- 
port at ion should be characterized by a suitable balance between the entanglement content and the degree of affinity 
to squeezed vacuum states. 

II. CV TELEPORTATION WITH SQUEEZED CAT-LIKE RESOURCES 

By replacing in Eq. ([3]) the two-photon term 1 1 , 1) 12 with the two-mode symmetric coherent state 17,7)12, with 
7 complex amplitude ((7,7)12 = |7)i <8> (7)2), the squeezed Bcll-likc state (j3|) is converted in the following squeezed 
cat-like superposition: 

|V)sc =AA^i 2 (C){cos ( 5|0,0)i2 + e je sin ( 5|7,7)i2}, (16) 

where the normalization factor is M = {1 + e - ' 7 ' sin 26 cos #} -1 / 2 . We remark that the possibility of engineering 
states of the form \ip)sc is confirmed by the recent successful experimental realization of their single- mode equivalent 
[23I ]. The characteristic function associated with Eq. (fl6|) . Xsc( a i, a 2) — sc(i J \Di(ai)D2(a2)\ip)sc, can be easily 
calculated, yielding a simple algebraic sum of Gaussian functions. We can thus compute analytically the fidelity of 
teleportation using squeezed cat-like states <|16|) as entangled resource, limiting the analysis, for simplicity, to the 
teleportation of input coherent states |/3). The maximization procedure is discussed in Appendix [A"l including the 
analytical expression for the fidelity. At fixed squeezing r, the optimal fidelity can be obtained by maximizing on the 
real amplitude I7I, as the only free parameter 

F op t{r) = maxjFsc (r, \~f\) . (17) 

It] 

In Fig. [4J we show the behavior of T ov t as a function of r for the input states |/3), with an optimized squeezed 
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FIG. 4: Optimal fidelity of teleportation T op t, as a function of the squeezing parameter r, with = 0. The fidelity corresponds to 
the teleportation of input coherent states \f3) with a squeezed cat-like state as entangled resource (dotted line). For comparison, 
we plot the fidelities obtained with a squeezed vacuum resource(full line) and with a squeezed Bell-like resource (dashed line). 
The inset gives the maximal value of the parameter 7I = 7 as a function of r. 

cat-like state as entangled resource. For comparison, also the optimal fidelities, corresponding to the teleportation 
with Gaussian squeezed state resources and with squeezed Bell-like states, are plotted. Optimized squeezed cat-like 
non-Gaussian resources yield a significant improvement of the fidelity with respect to the two-mode squeezed vacuum, 
but are less efficient than squeezed Bcll-likc resources. This result is in agreement with the behaviors of the von 
Neumann entropy, the non-Gaussianity, and the squeezed-vacuum affinity, reported in Figs. [5l and [5] panel I and 
panel II, respectively. In fact, the optimized squeezed cat-like superposition is less entangled, less non-Gaussian 
(according to the particular measure d n Q) 1 and less affine to the squeezed- vacuum, than the optimized squeezed Bell- 
like state. In particular, the optimized squeezed cat-like resource exhibits a marked reduction of Q, if compared with 
the corresponding plots shown in Fig. [6] panel II; in fact, Q keeps lower than 71%, for any r. 

III. CV TELEPORTATION WITH NOISY NON-GAUSSIAN STATES 

In this Section, we extend the procedures described above to the instance in which, while keeping the protocol 
ideal, the entangled resources are affected by thermal noise. In order to describe how a mixed resource arises from a 
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FIG. 5: Behavior of the von Neumann entropy E v n as a function of r, with |7| = 7 (see the inset in Fig. 2} for the optimized 
squeezed cat-like state (dotted line). For comparison, we also plot the quantities E v n corresponding to a squeezed vacuum (full 
line), and to a squeezed Bell- like state (dashed line). 
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FIG. 6: Behavior of the non-Gaussianity d„G (panel I) and of the squeezed-vacuum affinity Q (panel II) as a function of r, with 
I7I = 7 (see the inset in Fig. [4} for the optimized squeezed cat-like state (dotted line). For comparison, in panel I we plot d n G 
for the optimized squeezed Bell-like state (dashed line). 

pure one when a thermal background is turned on, we consider fields corresponding to independent contributions and 
converging on a common spatial volume. We follow a very effective scheme proposed by Glauber, based on the P- 
function formalism [25| . The superposition of two independently generated fields, both admitting a P representation, 
is simply the convolution of the individual P- functions: 



p(0 = J J d^iMe-O^CO- 



(18) 



As demonstrated in Ref. [33[ (see also Ref. [34[ ) , if at least one of the independent fields admits a valid P representation, 
then the normally ordered quantum characteristic function x_/v(a) of the superposition is given by the multiplication 
of the individual characteristic functions: 



Xn{o) = x { n '(«) X ( n («) • 



(19) 



Thus, the relation (|19p can be applied in the case of the superposition of a generic field p with a single-mode thermal 
field 



o {1) 



Hh 



(l + n th ) k + 



T \k)(k\ 



(20) 



whose normally ordered characteristic function is a Gaussian function of the form 



(21) 
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where n t h denotes the average number of thermal photons. Generalizing the formula (|19|) to the two-mode instance, 
we consider the superposition of a two- mode non-Gaussian state p = \ip)nG noi^l an d a two- mode thermal state p th = 

Pthi ® Pth2- Next, we specialize the analysis to the two- mode squeezed Bell-like state and to the two- mode squeezed 
cat-like state, given by Eqs. ([3]) and (fl6|) . respectively. The resulting field is then described by the symmetrically 
ordered characteristic function 

xi&Woa) = z~ nthMa ^~ nthMa * ? Xn G {ui,a 2 ) , (22) 

where n t h.i is the thermal parameter associated with mode i (i = 1,2), and XnG{pt\ 1 a 2 ) denotes the characteristic 
function associated with the considered non-Gaussian state. The state described by Eq. (f2"2"| is mixed as soon as 
n th,i 7^ 0. Of course, for 5 = we obtain the characteristic function for the superposition of a two-mode squeezed 

(2) 

vacuum state p = \Q (£| and a two-mode thermal state p th . As in the previous Sections, we can now exploit the 
mixed non-Gaussian resource, described by the characteristic function (j22|) . and use Eq. (|10|) to calculate the fidelity 
of teleportation, for instance in the case of input coherent states. The explicit analytical expressions for the fidelities 
of teleportation of input coherent states either with mixed squeezed Bell-like resources or with mixed squeezed cat-like 
resources, i.e. J-gg\r, n t h,i, nth,2, $) and J^s^ir, n t h,i, nth,2, M), are reported in Appendix [Al In the case of mixed 
squeezed Bell-like states, at given thermal numbers nth,i, nth,2, and fixed squeezing parameter r, the optimal fidelity 
is defined as: 

Fo P t(r,n thA ,n th ^) = max Tgg (r,n thA ,n th ^ 2 , 5) , (23) 



and the optimal angle 5 = 8 max is easily computed: 

1 / e~ 2r \ 

Smax = 7: arctan 1 + — ■ . (24) 

2 V l + nth, i + nth^' 

In the case of mixed squeezed cat-like states, at given thermal numbers nth,i> n th,2, and fixed squeezing parameter r, 
the optimal fidelity is defined as: 

T ov t{r,nth,\,nth,2) = max Tg^ \r, nth,u nth,2, M) ■ (25) 

In Eq. (|25p the maximization can be performed numerically. In Fig. we plot the optimal fidelities for the mixed 
non-Gaussian resources as a function of squeezing for various choices of the thermal parameters n t h,i = n t h,2 = nth- 
We observe that, as expected, the fidelity decreases for increasing nth'- the thermal noise sensibly reduces the success 
probability of teleportation. However, at fixed squeezing and thermal parameters, the non-Gaussian resources always 
perform better than the Gaussian twin beam. Moreover, also in the mixed instance, mixed squeezed Bell-like states 
perform better than mixed squeezed cat-like states. Importantly, for the chosen realistic values of n t h, the fidelity 
associated with the mixed non-Gaussian resources never drops below the threshold of classical teleportation with 
maximal fidelity T™^ x = 0.5. On the other hand, one can look for the threshold value n^ s \r) such that, at fixed 

r, one has T opt = F^Ls* ~ \- ^ n Fig- IH1 we plot n|^ s ' ) as a function of r. We sec that the threshold value nf h ls ^ 
associated to non-Gaussian resources is larger than that associated to Gaussian ones, and in any case sensibly larger 
than the reasonable realistic values we have considered. On the contrary, the fidelity associated with the mixed 
squeezed vacuum state, at low values of r, falls below the classical threshold, unless optimization is performed on the 
local degrees of squeezing, without modifying the Braunstein-Kimble protocol [35|. 



IV. CONCLUSIONS 



We have analyzed the performance of the Braunstein-Kimble continuous-variable quantum teleportation protocol 
using two particular classes of non-Gaussian resources, i.e. the squeezed truncated twin beam and the squeezed 
cat-like states. The extra freely tunable parameters admitted by these resources in comparison to standard Gaussian 
resources allow for an optimization procedure in order to maximize the fidelity of teleportation. For a given set 
of input states, such a procedure yields a sculptured output telcported state. In order to interpret the results and 
explain the enhancement of the fidelity, we have characterized and compared resources according to their degrees 
of entanglement, non-Gaussianity, and squeezed-vacuum affinity. We have shown that the optimized non-Gaussian 
resources possess both a high entanglement content and a high degree of squeezed- vacuum affinity. Finally, we have 
studied the efficiency of CV teleportation of input coherent states using, as entangled resources, mixed non-Gaussian 
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FIG. 7: Optimal fidelity T op t, associated with the teleportation of input coherent states \0) with mixed squeezed Bell-like 
resources (panel I) and mixed squeezed cat-like resources (panel II), as a function of the squeezing r, for several choice of the 
thermal parameters n t h,i = n t h,2 ~ nth- The curves representing the fidelities associated with mixed non-Gaussian entangled 
resources are plotted with the following plot style: nth = (full black line), n t h = 0.05 (long-dashed line), n t h = 0.10 (double- 
dot dashed line), and nth = 0.15 (dotted line). For comparison, in both panels, we also plot the fidelities associated to the 
mixed squeezed vacuum entangled resources, with n t h = (full gray line), nth = 0.05 (dashed line), nth — 0.10 (dot-dashed 
line), and n t h = 0.15 (long-dotted line). The inset in panel II gives the maximal value of the parameter | — y | = 7 as a function 
of r, with nth = 0, 0.05, 0.1, 0.15 (same plot style as for the fidelities). 




FIG. 8: Behavior of the threshold value nf h a as a function of r, for mixed squeezed Bell-like resources (dashed line), mixed 
squeezed cat-like resources (dotted line), and mixed squeezed vacuum (full line). 

states of the radiation field. We have shown that, although the thermal noise strongly affects the success probability 
of teleportation, the resources provided by the optimized mixed-state extension of pure non-Gaussian resources guar- 
antee sufficiently high fidelities, for realistic values of the average thermal photon number n t h and of the squeezing 
r. In particular, they assure better performances with respect to the mixed extensions of pure Gaussian twin-beam 
states. 

The present analysis can be generalized by considering more general non-Gaussian resources and fully realistic sit- 
uations that involve simultaneously nonidcal Bell measurements and dccoherencc induced by propagation in noisy 
channels. A further goal will be the study of optimal teleportation with non-Gaussian resources of two-mode and 
multimode states (36j . and the optimization of teleportation protocols with non-Gaussian resources with respect to 
local parameters and properties, extending the existing schemes for the local optimization of Gaussian resources, that 
yield a one-to-one correspondence between the content of entanglement present in a resource and maximal nonclassical 
output fidelity (3oT |. 
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APPENDIX A: ANALYTICAL EXPRESSIONS FOR THE FIDELITIES 



Here we report the analytical expressions of the fidelities of teleportation of input coherent states with the pure 
and the mixed non-Gaussian squeezed Bell-like states as resources (we do not report the whole calculation that is 
straightforward but somewhat cumbersome). The fidelity TsB{ r ,5) for the squeezed Bell-like state ([3]) reads: 



1 + e 2r + e Ar + e 2r cos 28 + (1 + e 2r ) sin 26 
^SB(r, S) = c - 2r(1 + e2r)3 " • ( A1 ) 



For S = 0, the state \ip)sB reduces to the twin beam, and Eq. (|A1[) gives the well known formula: 

T TwB (r) = T SB (r, 0) = 1 _ 2r . ( A2) 

1 + e zr 

The fidelity T { s%>{r, n t h, i, n t h, 2 f° r the mixed squeezed Bell-like state reads: 

rith), « 1 + e 2r f th + e ir f 2 h + e 2r f th cos 25 + [1 + e 2r f th ] sin 2,5 

where we have defined 

fth = 1 + n t h,i + n t h,2 ■ (A4) 

Of course, it is easy to check that, in the absence of thermal fields, i.e. n t h.i = n t h,2 — 0, we get fth = 1, and the 
fidelity Tgg reduces to Tsb- Notice that, in the limit of large thermal parameter, the fidelity of teleportation (|A3[) 
is strongly suppressed and vanishes asymptotically. 

Next we give the analytical expressions of the fidelities of teleportation of input coherent states with pure and mixed 
non-Gaussian squeezed cat-like states as resources. By putting, as usual, the phases 9 and </> to zero, the fidelity for 
the squeezed cat-like state (fll))) reads: 

cos 6 + e sin 6 + e 171 e 1 += 2r + e 1 +'' Jr sin S cos S 



f'sc( r i 5 il) = 7j T _ 2rVl T _ h i2 ■ ox , ■ (A5) 

(1 + e zr )(l + e n I sin 2d) 



It is worth noticing that, for 5^0 and/or 7^0, the fidelity (|A5j) reduces to the well known expression Ttwb(t), 
holding for a twin beam resource. A preliminary numerical optimization procedure for Eq. (|A5|) allows to fix the 
parameters arg 7 and S to the values arg 7 = and S — j, leading to a simplification of the fidelity to the formula: 



^o(r,| 7 |)= (i + e-t^T+l-hP) - ( A6 > 

Thus, at fixed squeezing r, the maximization can be carried out on the real amplitude I7I, as the only free parameter, 
i.e. J- op t = max| 7 | Tsc {f, \"f\)- In the particular instance of zero squeezing r — 0, it is easy to find value of the optimal 
fidelity T opt = [4(\/2 - l)]" 1 ~ 0.6035, for | 7 | = ln 1/2 (\/2 - 1)~ 2 ~ 1.3276. Considering the mixed instance, at fixed 

angle S — \ and phases <j) — ^, 8 = 0, arg7 = 0, the fidelity Jgg (r, n t h,i, nth,2, \ j\) for the mixed squeezed cat-like 
state reads: 



1 + e-^e^+^hh 

'-sc'{ r ' n th,i,n t h,2, = — 

where fth is defined in Eq. (|A4I 



T sc (r,n thA ,n th ,2,\l\) = e _ 2r(1 + + e _ W) , (AT) 
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